Magnetism and magnetotransport in disordered graphene 
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We perform Monte Carlo simulations to study the interplay of structural and magnetic order in 
single layer graphene covered with magnetic adatoms. We propose that the presence of ripples in 
the graphene structure can lead to clustering of the adatoms and to a variety of magnetic states such 
as super-paramagnetism, antiferromagnetism, ferromagnetism and spin glass behavior. We derive 
the magnetization hysteresis and also the magnetoresistance curves in the variable range hopping 
regime, which can provide experimental signatures for ripple induced clustering and magnetism. 
We propose that the magnetic states in graphene can be controlled by gate voltage and coverage 
fraction. 

PACS numbers: 73.20-r,73.20.Hb,75.20.Hr 



I. INTRODUCTION 

Graphene is probably one of the most remarkable dis- 
coveries in condensed matter physics in the last decadei^i^. 
The material is a two-dimensional (2D) crystal composed 
of carbon (C) atoms with sp^ hybridization, that is, 
graphene is one atom thick and thus the thinnest cloth 
in nature. Because of its low dimensionality, it does not 
show structural long-range order in its free form, as per 
the Hohenberg-Mermin- Wagner (HMW) theorem^^, but 
it can present a flat phase at low temperatures due to 
non-linear effects or in the presence of a substrate, scaf- 
folds, contacts, or impurities that break explicitly the 
translational symmetry perpendicular to the graphene 
planed. The leftovers of the fluctuations that forbid long- 
range order are found in the form of frozen ripples in 
suspended^'^ samples. 

While it was theoretically predicted early on that 
graphene by itself would not be magnetic^, it has been 
shown that adatoms in the graphene surface can eas- 
ily form magnetic moments due to graphene's unusual 
electronic properties such as low density of states and 
chirality-. Moreover, because graphene has a low den- 
sity of states close to the Dirac point, the Kondo effect 
is suppressed^iiS, allowing for the appearance of mag- 
netic states. On the other hand, because of its low di- 
mensionality, long-range magnetic order is inhibited, and 
the intrinsic coupling between structural and magnetic 
fluctuations may lead to very inhomogeneous spin tex- 
tures (such as of Griffiths phasesii) that can be found in 
complex itinerant magnetic systems such as disordered 
Kondo Lattices^. 

In this paper, performing a series of Monte Carlo simu- 
lations over disordered realizations of magnetic adatoms 
in a graphene sheet, we examine the magnetic correla- 
tions of the local moments in the magnetization and in 
the magnetoresistance curves. In the case of adatoms 
such as Hydrogen (H), where the probability of adsorp- 
tion changes substantially according to the local curva- 



ture of the graphene sheet, we show that the formation 
of local magnetic moments leads to an interesting inter- 
play between the correlation due to the RKKY interac- 
tion and the ripples, generating a variety of magnetic 
textures, tendency to clustering at low concentrations 
and a percolative transition at higher concentration of 
adatoms. Since the overall magnetic response of a single 
layer is rather small compared to the usual response in 
bulk, magnetotransport measurements are probably the 
easiest way to probe magnetic correlations in graphene. 
Based in our Monte Carlo results for the magnetization, 
we calculate the magnetoresistance curves in graphene 
for strong disorder. These curves can offer clear experi- 
mental signatures for the presence of macroscopic mag- 
netic states in graphene in the regime of variable range 
hopping^'^. 

The structure of the paper is as follows: in sec. II 
we introduce the spin Hamiltonian for the magnetic 
adatoms; in sec. Ill we derive the heuristic rules for 
adsorption of adatoms in graphene, comparing different 
realizations of disorder for different adsorption probabil- 
ity distributions, which we define in terms of the cur- 
vature (height) of the ripples. In sec. IV, we derive 
the magnetization curves and in sec. V we calculate the 
magnetoresistance. Finally, in sec. VI we present our 
conclusions. 



II. RKKY HAMILTONIAN 

Our starting point is the tight-binding Hamiltonian of 
the electrons in graphene (we set H = 1)^: 

Htb = -t5]5][4(R.K(Rj)+h.c.], (1) 

ihj) 

where ao-(Ri) {b^(Ri)) annihilates and electron with 
spin a —^,1 on sublattice A (B) at position R^, and 
t (« 2.7 eV) is the nearest neighbor hopping energy. 
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When adatoms are added to graphene they can hy- 
bridize with an energy V to the C atoms and if the local 
Coulomb energy U is sufficiently large, a local moment 
of spin S is formed^ at a site R^. This spin Si interacts 
with the graphene electrons via an exchange interaction, 
Jk Ri —V'^/U, which is described by the Hamiltonian: 

= Jfe^^,-s,, (2) 

i 

where Si is the graphene electron spin. Eq. ([1]) together 
with ^ describe a Kondo lattice in graphene. 

The Kondo interaction ^ induces an indirect kinetic 
exchange interaction between local moments, the RKKY 
interaction, which depends on the chemical potential fi 
and hence can be controlled with an external gate volt- 
age. Therefore, the nature of the magnetic states in 
graphene can be controlled by the application of a trans- 
verse electric field, a situation that never occurs in met- 
als. For low carrier concentrations, i.e. close to the Dirac 
point, the interaction between spins located in the same 
sublattice is ferromagnetic, xaaCR) = Xsb(R.) oc 
while it is antiferromagnetic if they belong to opposite 
sublattices (xab(R) oc -l/R^)l±d^J^. Given the po- 
sition of the magnetic moments and the effective interac- 
tion between them, we can perform numerical simulations 
to obtain the magnetic properties of the system. 

When an adatom sits on top of a carbon atom, the 
carbon sp^bonds are locally distorted and acquire a sp'^ 
character as in diamond, with the carbon underneath 
the impurity being pulled out of the plane^^. Although 
a flat graphene sheet has almost no magnetic anisotropy 
due to the very small spin-orbit couplingiS., in the pres- 
ence of ripples?", lattice distortions of the sp^ type that 
are generated by the presence of adatoms such as H can 
substantially increase the spin-orbit interaction close to 
the adatoms, generating a local spin-orbit coupling up to 
Ago « 7 meV— , which corresponds to an out of plane 
magnetic anisotropy of « 7 T. In the range of tempera- 
ture ksT < A50, the minimal model for the interaction 
between adatoms in graphene should be correctly cap- 
tured by the physics of the Ising model: 

HeS = J2 JKKKY{nj)S^S^ - g^lBH,., ^ . (3) 

where Jrkky {fij ) is a spatially dependent RKKY inter- 
action that depends on and Hf,x is an external magnetic 
field. For fi ^ t, the period of the RKKY oscillation is 
long compared to the typical atomic distances and we 
only need to consider the power-law decay of the inter- 
action. The exchange is given by JAAi^) = JBsir) = 
Jo exp{—rij / rg) / r^j for spins in the same sublattice and 
JAsir) = —Jq eTq){—rij/ro)/rfj for spins in different sub- 
lattices, rij is given in units of the lattice constant 
a — 2.46 A and tq w 12 A is a cut-off distance that 
is introduced to limit the range of the parameters in the 
numerical calculations. 



III. RIPPLES AND ADSORPTION 
PROBABILITY 

In a flat graphene layer, the magnetic interaction be- 
tween impurity spins would generate uniform magnetic 
states. Nevertheless, ripples in the graphene structure 
break the translational symmetry and lead to a inhomo- 
geneous situation where the adatoms have preferential 
sites for hybridization. Due to ripples, the topography 
of the surface in graphene is not flat but has curvature, 
which appears in in the form of "valleys" and "hills" . In 
graphene, the typical ratio between the height (h) and 
the lateral size (L) of the ripples is h/L « 0.1 — 0.2^. On 
top of a hill, the characteristic curvature of the ripples 
distorts the sp^ bonds to an extent that they acquire sp^ 
character, with a significant hybridization between the a 
and TT bands. In this situation, the adsorption of a H 
atom on a hill-C has been shown to be energetically fa- 
vorable, helping to stabilize the ripple^^. On the other 
hand, for perfectly flat graphene, the adsorption of a H 
atom costs some energy to locally distort the sp'^ bonds 
and pull the hydrogenated C atom out of the plane. The 
difference in the chemisorption energy between the two 
cases, recently calculated by ab initio methods, can be 
as large as 2.5 eV. These results also indicate the exis- 
tence of a minimal curvature for the ripples {h/L > 0.12) 
above which the adsorption of a H atom can be strongly 
favorable^-. The substantial change in the chemisorp- 
tion energy of H with the local curvature in graphene is 
qualitatively consistent with previous theoretical — and 
experimental^i results in nanotubes, where the binding 
energy has been shown to change dramatically with the 
radius of the tube. 

Assuming that graphene is supported on a substrate, 
only one surface is available for hybridization. Notice 
that because of curvature, the Pz orbitals that arc lo- 
cally perpendicular to the graphene surface, approach 
each other in the valleys but distance themselves in the 
hills. Contrary to the hill-C case, the local curvature 
is expected to inhibit the H adsorption in the valleys, 
since in order to stabilize the C-H bond relative to the 
other three C sp'^-like bonds, the H has to pull the C 
atom out of the plane against the curvature, at the ex- 
pense of an additional cost of energy, compared to the flat 
case (the difference in chemisorption energies between the 
hill-C and the valley- C cases could easily amount to sev- 
eral eV). The H adatoms therefore hybridize much more 
easily in the hills than in the valleys, leading to a "per- 
colative" structure as the planar density of adatoms is 
increased. Therefore, the probability that an adatom hy- 
bridizes with a hill-C is larger than a valley-C atom. 

Let us assume, for the sake of the argument, that 
the ripples in graphene have a simple sinusoidal form 
(see ref. HI), that is, if z is the perpendicular coordi- 
nate of the C atom relative to the flat situation, then 
z{x, y) = A sin{kxX + kyy -\- 4>) where A, (j), kx and ky are 
random parameters. A large number (up to 200) of such 
sinusoidal waves can be superimposed to obtain a ran- 
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FIG. 1: (a) Graphene sheet with random ripples. The color 
map represents the height increasing from blue to red in the 
interval z G [— Zmax, Zmax]- (b) -(c) graphene sheet covered 
with adatoms for a lower cut-off fto = 0.28 Zmax and /lo = re- 
spectively (see main text), (d) Zoom- in showing the adatoms 
on top of a hill (red rectangle). 



domly curved sheet, as shown in Fig. [T] (a) for a system 
with 2 X 600 X 600 C atoms. Next, we decide a strategy 
to allow the incorporation of adatoms by the graphene 
sheet for some quenched realization of ripples. 

Given the existence of a lower bounded range of lo- 
cal curvature where the adsorption is most likely to 
happenS^, and also the fact that the adsorption is much 
easier on the top of a hill rather than in the bottom of a 
valley in graphene, we assume a probability distribution 
that grows monotonically with the height of the ripples 
and has some effective lower cut-off, /iq: below which the 
adsorption is unlikely to happen. As a toy model, we 
may assume for instance that the probability of adsorp- 
tion varies linearly with z through the following criterion: 
we allow the incorporation if 



z > ho + hi 



(4) 



where hi is a random number varying from to Zmax^h-Q, 
and ho G [—z^ax, Zmax] is the lower cut-off of the distri- 
bution, with Zmax as the maximum height. The choice 
of another distribution function such as, for instance, a 
step function, 0{z — /iq), or some other distribution with 
a finite tail for z < h^, should lead to similar conclusions 
regarding the physical properties, as far as the probabil- 
ity of adsorption on the top of the hills is large compared 
to the probability adsorption in the valleys. 

If the number of atoms which are available for adsorp- 
tion is fixed, by varying ho we obtain different incorpora- 
tion fractions, as shown in Figures [1] (b),(c). Notice that 
for a given realization of ripples, as one increases the 
adatom concentration, one obtains a percolative struc- 
ture. At low coverage densities the adatoms form clus- 
ters. Hence, ripples naturally lead to clustering. We 
can also fix the total number of adatoms attached to the 
sample and see how the different probability functions 
change the coverage structure of the system. In this case, 
the overall structure is the same of Figures [T] (b),(c) but 



different probability distributions lead to different con- 
centration of adatoms on top of the hills, as it will be 
clear in the next section. 



IV. MAGNETIC TEXTURES 

The details of the clustering depend on the nature of 
the adatoms used. For instance, if one uses H atoms, one 
might imagine that only H atoms that are separated by 
more than one lattice spacing are stable since H atoms 
which are in neighboring C sites can recombine into II2 
molecules and leave the graphene surface^. Hence, we 
have considered two different situations, namely, either 
the adatoms have nearest neighbors or not. We studied 
the location of magnetic adatoms using the same incorpo- 
ration rules discussed above, with and without allowing 
the existence of nearest neighbors. For the system with- 
out first neighbors, we randomly remove the neighbors, 
in order to avoid any artificial distribution of spins. We 
found that the magnetic response of these two systems 
are rather different and can be distinguished experimen- 
tally. 

For a graphene sheet with full coverage, as shown in 
Fig. [21(a), we find that the moments order antiferromag- 
netically, as expected. Nevertheless, when we associate 
the local probability of adsorption of the adatoms to the 
ripples and we start to dilute the coverage for increasing 
values of /iq, the situation is rather more complicated. 
What we observe in figure [5] (b)-(d) is the destruction of 
the antiferromagnetic order on the clusters that are lo- 
cated on top of the ripples. The system becomes highly 
frustrated and finally for very low coverage we find that 
the system consists basically of weakly interacting iso- 
lated moments and clusters of few spins. One notices 
that the destruction of antiferromagnetic order is accom- 
panied by a ferromagnetic tendency, that is, the clusters 
tend to either remain antiferromagnetically ordered and 
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FIG. 2: Spin configuration on top of a hill for (a) homogeneous 
coverage and (b) -(d) ho — -1, and 0.28 Zmax with adatom 
concentrations x = 0.5, 0.11 and 0.02 respectively. In black, 
spin up; in orange, spin down. (fcsT = O.OlJo). 
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FIG. 3: Left: (a) Magnetization and (b) staggered magnetiza- 
tion as a function of an external magnetic field for a fully cov- 
ered graphene sheet (black circles) and for ho—O (red squares) 
with ksT = O.OlJo. 



interact ferromagnetically with nearest clusters (see Fig. 
[5] (c) ) or order ferromagnetically in a situation that is 
reminiscent either of super-paramagnetism or a spin glass 
stated. "In this paper we have considered only the av- 
erage magnetic properties while disregarding rare events 
that can lead to Griffiths phases ??. However, rare event 
physics can play an important role in these disordered 
magnetic systems and may contribute to the magnetic 
response. This is a problem, however, beyond the scope 
of the current manuscript. 

For the configurations discussed above, we calculate 
the total average magnetization M = Ma + Mb and the 
average staggered magnetization, Ms — \Ma — Mb\, in 
terms of the magnetization of each sublattice. Ma and 
Mb , as a function of the magnetic field Hex ■ AH the nu- 
merical simulations were performed for an Ising Hamil- 
tonian with 2 x 600 x 600 sites. We used 10000 Monte 
Carlo steps for warm up and another 10000 steps to col- 
lect the data. The magnetization is computed as an av- 
erage over 160 to 400 realizations. As can be seen in Fig. 
[3]Ja) and (b) the maximally covered honeycomb lattice 
has always an antiferromagnetic correlation between the 
spins but as the coverage is reduced, this correlation is 
strongly suppressed, changing the shape of the hysteresis 
loops. 

If we suppress the existence of nearest neighbors, one 
favors ferromagnetic interactions, once the interaction 
between next-nearest neighbors is ferromagnetic. How- 
ever, as the magnetic atoms are randomly removed, we 
never produce a regular lattice of spins and do not have 
a perfect ferromagnetic system. Instead, for the sys- 
tem that was originally a regular honeycomb lattice, we 
see in figure [3] (a) the formation of magnetic domains 
in the form of antiferromagnetically correlated ferromag- 




FIG. 4: Spin configurations for the same hill considered in 
figure [2] but with the random removal of nearest neighbors, 
(a) homogeneous honeycomb lattice, and (b) -(d) ho = —1, 
and 0.28zmax respectively. The adatom concentrations are 
X — 0.24 in the homogeneous lattice (a,), x = 0.22 (b), x = 
0.08 (c) and x — 0.01 (d). In black, spin up; in orange, spin 
down (ksT = O.OlJo). 



netic stripes. In the dilute situation [see Fig. |4l^b)-(d)] 
these ferromagnetic stripes are broken apart leaving be- 
hind ferromagnetic strings. The situation is clearly quite 
complex from the magnetic point of view and it is im- 
portant to distinguish these various spin textures exper- 
imentally. 

In the case of an originally regular honeycomb lattice 
without neighboring adatoms, the magnetization M{H) 
in Fig. [5] shows a large coercivity due to the ferromag- 
netic stripes (black circles) , which decreases with smaller 
incorporation fractions. Whereas for large incorporation 
fractions we still can see the signature of antiferromag- 
netic correlations in the shape of the hysteresis curve, for 
smaller incorporation fractions the systems consists of 
ferromagnetic clusters with a typical ferromagnetic hys- 
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FIG. 5: Magnetization as a function of magnetic field for 
a lattice without first neighbors: homogeneous lattice (black 
circles), ho=0 and fcsT = O.OlJo (red squares) and ksT — 
0.05 Jo (green triangles). Inset: ho=0 and ksT = 0.2 Jo (blue 
diamonds) . 





FIG. 6: The three columns show the spin configuration (top), the details of the configuration for a small region of the system 
(middle) and the total magnetization versus magnetic field for (a) ho — 0.28 Zmax (b), ho — and (c) a random distribution. 
The three systems contain the same concentration of adatoms, x = 0.05. 



teresis loop (green triangles and red squares). By increas- 
ing temperature, we observe a decrease in the coercivity 
of the hysteresis loop. For ksT > 0.2Jo (see inset of 
Fig. [5]), the spins interact weakly and we obtain a hys- 
teresis loop that resembles a paramagnetic behavior. In 
this case, for low magnetic fields, H^x, we find a universal 
hnear dependence in the magnetization, M oc Hex- 
It is important to point out the differences between the 
type of magnetic dilution we present here, compared to 
the usually magnetic diluted lattice, where the spins are 
randomly located, without any preferential position. In 
honeycomb lattices, randomly diluted spin systems have 
a site percolation transition at pc = 0.6970422.. For near- 
est neighbor ferromagnetic interactions, the occurrence 
of a magnetic transition as a function of the dilution p 
has been shown to coincide with the percolation transi- 
tion p^. For the antiferromagnetic case, recent results 
for the quantum Heisenberg Hamiltonian with site dilu- 
tions^ found that the magnetic long range order persists 
above Pc- In any case, in the randomly diluted system 
the moments are always isolated and weakly coupled at 
very small concentrations {p <C Pc)- In this work, we 
have shown that due to the distribution of spins accord- 
ing to the structure of the ripples, we suggest the possible 
existence of strongly coupled magnetic clusters on top of 
the highest graphene hills, even at very low concentra- 
tion of adatoms. Also, due to the ferromagnetic coupling 
between next nearest neighbors, the moment of a given 
cluster will be bigger than zero even for a mostly antifer- 
romagnetic one. 

In Figure [SI we compare the coverage structure and 



magnetization curves for three systems with the same 
adatom concentration, x = 0.05, but different incorpora- 
tion probability distributions. The total concentration 
is the same in the three systems but ho sets regions 
where the presence of adatoms is forbidden. As a result, 
we can see the presence of strongly correlated clusters 
for ho — 0.28 Zmax, while for a random distribution of 
adatoms we found almost isolated spins. In addition, one 
can see the clear difference between the hysteresis loop 
in the correlated case (ripples), shown in Fig. [6] (a), (b), 
compared to the completely random one [Fig. [S] (c)] . In 
the correlated case, we found regions with high concen- 
tration of spins and very strong antiferromagnetic corre- 
lations between them at low temperatures (fcsT <C Jo), 
in contrast with the random case, where both hysteresis 
and the coercitive fields are very small. 

V. MAGNETORESISTANCE 

The interplay between the corrugated nature of 
graphene and its magnetic properties, which gives rise 
to the spin textures discussed above, can be probed 
by magneto-transport measurements. It was recently 
demonstrated that the incorporation of adatoms can dra- 
matically change the transport properties in graphene^, 
where the system can go from a metallic behavior to a 
variable range hopping regime. In the case of magnetic 
impurities, due to the exchange interaction between the 
spin of the carriers and the adatom localized moments, 
the transport properties will be further modified by the 
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interaction between spin and charge degrees of freedom. 
For strongly disordered graphene, which characterizes 
the regime of variable range hopping, the carriers are 
trapped with a binding energy Ei. Even at Hex — 0, the 
exchange interaction fully polarizes the localized spins in- 
teracting with a given localized carrier and increases its 
binding energy^''. 

If one considers the localized spins in the mean-field 
approximation, the exchange energy is given by Eg^ = 
JkMg(s) where is the exchange interaction between 
localized spins and carriers, is the saturation mag- 
netization of the magnetic impurity spins and (s) is the 
mean value of the localized carriers spin. The total bind- 
ing energy in the presence of the exchange interaction is 
given by Eb = Ei + Eex and has an associated Bohr ra- 
dius ^. It this situation, each spin is fully polarized in a 
random direction. If we apply an external field i/ex, the 
localized spins in the whole sample begin to align with 
the field and because of the exchange interaction between 
them and the band carriers, H^x gives rise to a splitting 
of the conduction band. As a consequence, there is a de- 
crease in the binding energy by AE^x = JkM{T^ H){s). 
As Eh oc the magnetic field produces an increase in 
the effective Bohr radius 

CesiT,H) = {E,+Jk{s)M,[l-\M{T,H)\/M,])-'/\ (5) 

This effect is similar to the one observed in magnetically 
semiconductors as EuS or CdMnTe^. 

The hopping probability between two states at a dis- 
tance r is then given by P = exp(— 2r/^cff — W/ksT) 
where W is the energy difference between the two states 
and ^eff is the characteristic Bohr radius of the localized 
states under the effect of an external magnetic field. Fol- 
lowing the original Mott derivation, the resistance is^, 

p = Poexp(To(ff)/r)i/3, (6) 

where 

To{H) cx 13.8/ {kBN{fi)e,s) , (7) 

and Af(/i) is the density of states at the Fermi energy. 
Since Tq (x 1 — \M{T, H)\/Ms for small magnetization, we 
can extract the magnetoresistance curves of the system 
from our previous Monte Carlo results. In all curves, we 
normalize the scale of the magnetoresistance by the zero 
field case, R{Hex = 0) = 1. From the experimental point 
of view, the change in the temperature and gate voltage 
parameters in the range where Eq. ([5]) is applicable will 
additionally rescale the magnetoresistance curves. The 
scaling analysis of the magnetoresistance provides addi- 
tional information allowing the empirical determination 
of the exchange interaction Jk and the binding energy 
E,. 

We begin our analysis by considering non-interacting 
spins. In this case, the response is paramagnetic and the 
magnetization curves follow a typical Brillouin function 

2S + 1 /'2S+1 \ 1 / 1 \ 
MiHex,T) oc ^^coth(^^^.j _ _coth(^-.j 

(8) 



where x = g^BSHex/kBT, and S is the spin quan- 
tum number. Assuming g = 2 and S — 1/2, we show 
in Fig El^a) the magnetoresistance curves for a system 
with non-interacting spins at different temperatures. The 
magnetoresistance is negative, reaching a minimum value 
when all the spins are aligned by the magnetic field. For 
a fixed magnetic field, the resistance grows with temper- 
ature. 

Next, we present the magnetoresistance curves corre- 
sponding to the simulations in the previous sections. If 
we allow the presence of neighboring magnetic adatoms. 
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FIG. 7: (a) Magnetoresistance for non-interacting spins at 
three different temperatures: fcsT = 0.01Jo(solid line), 
ksT = O.lJo (dashed line) and ksT = 0.2Jo (dotted line), (b) 
Magnetoresistance for the three different incorporation distri- 
butions shown in Fig|B] /io=0.28 Zmax (black circles), /io=0 
(red squares) and a random distribution (blue diamonds), (c) 
Magnetoresistance in the case without first neighbor H atoms 
for /io=0 and different temperatures: ksT = O.OlJo (red 
squares), fcsT = 0.05Jo (green triangles) and ksT — 0.2Jo 
(blue diamonds). All curves are normalized by the corre- 
sponding zero field resistance R{Hex = 0) = 1 and are sym- 
metric for Hsx < 0. 
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the antiferromagnetic correlations between neighbors is 
strong and the minimum value of the magnetoresistance 
is only obtained for very large values of the magnetic 
field, set in energy units of ksT. In Fig. [7fb) we com- 
pare the magnetoresistance curves that correspond to the 
configurations shown in Fig. [6] for a fixed concentration 
of adatoms, x = 0.05, and different incorporation proba- 
bihties. In the correlated case, [Fig.[nja),(b)], because of 
the fixed concentration of adatoms, the top of the highest 
ripples accumulates larger antiferromagnetic clusters for 
larger values of /iq, requiring the application of stronger 
magnetic fields to decrease the resistivity. On the other 
hand, for the case of random distribution of spins, the 
antiferromagnetic correlations are much weaker. Notice 
that for systems with hysteretic behavior in the magneti- 
zation, M is zero at the coercive field. Consequently, the 
magnetoresistance at the coercive field is higher than its 
value for H^x = 0, giving rise to a positive magnetoresis- 
tance at small fields. 

In the scenario without nearest H neighbors, we con- 
sider the probability distribution with /iq = 0, for dif- 
ferent temperatures [see Fig[71[c)]. For ksT > 0.2 Jq the 
hysteretic nature of the magnetoresistance is consider- 
ably reduced and we obtain a curve that is similar to the 
case of non-interacting spins. Neverthless, we note that 
the width of the resistivity peak is anomalously small, in- 
dicating that the effective size of such "non-interacting" 
moments is considerably larger than the moment of iso- 
lated adatoms. This is a clear indication for the presence 
of small magnetic clusters in graphene. For even higher 
temperatures, ksT ^ Jq, the resistance decreases with 
H following p oc exp(— ai/), where a is a constant that 
depends on the exchange interaction, density of states 
and temperature. For weak external fields H^x Ago, 
where A50 is the spin orbit coupling, the magnetoresis- 
tance will be strongly anisotropic and can vanish when 
the magnetic field is parallel to the graphene plane. 



VI. CONCLUSIONS 

In summary, we have discussed the possible magnetic 
states of magnetic adatoms on top of a rippled graphene 
sheet. We proposed the scenario where the non-trivial to- 
pography of graphene correlates the local moments of the 
adatoms, generating clusters and non-trivial macroscopic 
magnetic states. The magnetic order can be very sensi- 
tive to the value of the chemical potential, the adatom 
coverage, and the "repulsion" between adatoms for near- 
est neighbor atoms. 

We have found that while a perfectly covered graphene 
sheet has a strong tendency towards antiferromagnetism, 
the presence of ripples leads to more complex magnetic 
textures, with increased ferromagnetic tendency and even 
glassiness. We show that the hysteretic behavior of the 
magnetization provides a simple way to study these mag- 
netic orderings. Nevertheless, the presence of different 
magnetic structures affect directly the transport proper- 
ties. In the variable range hopping regime, the system 
presents a universal negative magnetoresistance that de- 
pending on the ratio between the exchange interaction, 
temperature and coverage fraction can have also show 
hysteretic behavior. The magnetotransport curves in 
graphene can provide clear experimental signatures for 
these non-trivial magnetic states. 
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